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Abstract 

We review recent results on higher-order calculations to squark and gluino production and decay at the LHC, as 
obtained within the Collaborative Research Centre / Transregio 9 “Computational Particle Physics“. In particular, 
we discuss inclusive cross sections, including the summation of threshold corrections, higher-order calculations for 
specific squark production channels and for top squark decays, and next-to-leading order calculations for exclusive 
observables matched to parton showers. 

Keywords: supersymmetry, higher-order QCD corrections 


1. Introduction 

After the discovery of a Higgs particle in the first 
phase of the LHC in 2012, the key scientific goal for 
LHC run 2 starting in 2015 is to search for and explore 
physics beyond the Standard Model (SM). Supersym¬ 
metric (SUSY) theories m are among the most attrac¬ 
tive extensions of the Standard Model. They allow for 
the unification of the electromagnetic, weak and strong 
gauge couplings 1301 and provide a solution to the fine- 
tuning problem of the SM Gna, if some of the super- 
symmetric particles have masses near the TeV scale Q 
Furthermore, in supersymmetric models with R-parity 
conservation (am, the lightest supersymmetric particle 
(LSP) is stable and may constitute the dark matter in the 
universe ifTOlfTTI . 

A generic LHC signature for supersymmetric mod¬ 
els with R-parity conservation are cascade decays of the 
strongly interacting supersymmetric particles, squarks 
and gluinos, which terminate in a weakly interacting 
LSP and thus result in missing energy signatures. Cur¬ 
rent LHC searches for supersymmetry in jets plus miss¬ 
ing energy final states place lower limits on the masses 
of squarks and gluinos at around 1.5 TeV Note, 

however, that these limits are not generic for supersym¬ 


metry, but depend on certain assumptions about the su¬ 
persymmetry breaking and the resulting SUSY particle 
mass spectrum. 

Given the importance of SUSY searches at the LHC, 
accurate theoretical predictions for the production and 
decay of supersymmetric particles are required. The 
precision calculations are crucial to derive accurate lim¬ 
its on SUSY masses and couplings, or to determine the 
properties of supersymmetric particles in the case of dis¬ 
covery. 

In this contribution we shall review the calculation 
of higher-order QCD corrections to squark and gluino 
production and decay at the LHC. In Sectionwe first 
present results for the inclusive production cross sec¬ 
tions, including the summation of large logarithmic cor¬ 
rections. The higher-order QCD corrections to SUSY 
particle decays are discussed in Section Fully differ¬ 
ential next-to-leading order predictions for squark pro¬ 
duction and decay, including the effects of parton show¬ 
ers, are presented in Section We conclude in Sec¬ 
tion |5] 
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2. Squark and gluino production at the LHC 

We consider the minimal supersymmetric extension 
of the Standard Model (MSSM) d [13 where, as 
a consequence of /^-parity conservation, squarks and 
gluinos are pair-produced in proton-proton collisions at 
the LHC: 

pp qq,q^,qg,ggjJiJ2h + X. ( 1 ) 


will thus focus on threshold resummation and the im¬ 
pact of non-degenerate MSSM particle spectra on the 
QCD corrections. The summation of the next-to-leading 
logarithmic corrections and the corresponding tool Nll- 
FAST are presented in Section |2.1[ In Section |2.2| we 
comment on the recent resummation of NNLL terms, 
while the NLO-QCD effects for generic MSSM spectra 
are discussed in Section [231 


The production of top squarks (stops), has to be 
treated separately, because the strong Yukawa coupling 
between top quarks, stops and Higgs fields gives rise to 
potentially large mixing effects and mass splitting US). 
In Eq. 0 and throughout the rest of this paper, t\ and 
t 2 denote the lighter and heavier stop mass eigenstate, 
respectively. For the other squarks we suppress the chi¬ 
ralities, i.e. q = (§L, §/?), and do not explicitly state the 
charge-conjugated processes. 

The cross sections for the hadro-production of 
squarks and gluinos in the MSSM are known includ¬ 
ing next-to-leading order (NLO) QCD |[T7H2T]| and elec- 
troweak ESHSIl corrections. The QCD corrections are 
particularly significant and included in the public com¬ 
puter code Prospino |[T7UT9l [32l. A large part of the 
QCD corrections can be attributed to the emission of 
soft gluons and can be taken into account to all orders 
in perturbation theory by means of threshold resum¬ 
mation techniques at next-to-leading logarithmic (NLL) 
or next-to-next-to-leading logarithmic (NNLL) accu¬ 
racy liSSHSOl . NLO-rNLL predictions for the production 
of strongly interacting MSSM particles as implemented 
in the computer code Nll-fast (331 123 (221 ESI SOI are 
currently state-of-the-art, and are employed by the LHC 
experiments and by large parts of the theory commu¬ 
nity to interpret search limits and constrain the MSSM 
parameter space Eiiisa. 

The calculations implemented in Prospino and Nll- 
fast assume that five flavors of left- and right-chiral 
squarks, d^R, Cl,r, Sl,r, and bi^R, are mass- 

degenerate. NLO-QCD predictions for generic MSSM 
spectra have been presented in the literature re¬ 
cently (53^5711 . and we shall comment on those in Sec¬ 
tion |2]3 Moreover, with the MadGolem EHEllEll 
and MadGraph5_aMC@NLO (601 frameworks there 
exist automated tools for the calculation of NLO-QCD 
corrections to generic supersymmetric processes at the 
LHC. 

While the NLO-QCD corrections for degenerate 
squarks have been available for many years (ItKISI, the 
summation of large logarithmic threshold corrections 
and the extension of the NLO calculations to generic 
MSSM spectra have been achieved more recently. We 


2.7. NLL threshold re summation 

A significant part of the NLO-QCD corrections to 
heavy particle production at the LHC can be attributed 
to the kinematic region where the partonic center-of- 
mass energy is close to the production threshold. Near 
threshold, the NLO corrections are typically large, with 
the most significant contributions coming from soft- 
gluon emission off the colored particles in the initial and 
final state. The contributions due to soft gluon emis¬ 
sion can be taken into account to all orders by means 
of threshold resummation. Threshold resummation for 
MSSM squark and gluino pair-production processes at 
NNL accuracy is described below, following the work 
presented in Refs. (33l|35l|371|38l|40l . We shall com¬ 
ment on NNLL resummation in Section [231 

The resummation for (2 ^ 2) QCD processes has 
been studied extensively in the literature, specifically 
for heavy-quark (^ [62l and jet production (6314651 . 
The calculations presented in Refs. (33l [^ [JTl [^ [40l 
make use of the framework developed there. 

The hadronic threshold for the inclusive production 
of two final-state particles k,l with masses and 
mi corresponds to a hadronic centre-of-mass energy 
squared that is equal to S = {m^ + m/)^. Thus we define 
the threshold variable p, measuring the distance from 
threshold in terms of energy fraction, as 

{nik + mif 

P = -?-• (2) 


The numerical results presented below are based on the 
following expression for the NLL-resummed cross sec¬ 
tion, matched to the exact NLO calculation (17] [TSl : 

NLO+NLL/^ fix ,2\ ^NLO / tlx ,2\ 
o-pp^ki },p ) = o-pp^j^iipAm },p ) 




dNp-^ fi/p(N + l,p^)fj/p(N + l,p^) 


X P^) Inlo ] • (3) 


The fnp.fjip are the parton distribution functions in 
Mellin space, and the last term in the square brack¬ 
ets denotes the NLL resummed expression expanded to 
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NLO. fi is the common factorization and renormaliza¬ 
tion scale. The resummation is performed after taking a 
Mellin transform (indicated by a tilde) of the cross sec¬ 
tion, 


d-pp^ki{N,{m^])= f dp ^ (Tpp^kiip,{m^\)- ( 4 ) 
Jo 

To evaluate the contour CT of the inverse Mellin 
transform in Eq. 0 the so-called “minimal prescrip¬ 
tion” IMIl has been adopted. The NLL resummed cross 
section in Eq. 0 reads 

1 

xAi(N + 1 , Q^,pAAj(N + 1 , Q\pA 

xdll^ki,i(N+hQ\pA, ( 5 ) 

where are the color-decomposed leading-order 

cross sections in Mellin-moment space, with I labelling 
the possible color structures. Here the hard scale = 
(rrik + mif has been introduced. The perturbative func¬ 
tions Cij^kij contain information about hard contribu¬ 
tions beyond leading order. This information is only rel¬ 
evant beyond NLL accuracy and therefore Cij^kij = 1 
is used in the calculations. The functions A/ and Aj sum 
the effects of the (soft-)collinear radiation from the in¬ 
coming partons. They are process-independent and do 
not depend on the color structures. They contain the 
leading logarithmic dependence, as well as part of the 
subleading logarithmic behaviour. The expressions for 
Ai and Aj can be found in the literature |[35]| . The resum¬ 
mation of the soft-gluon contributions, which depends 
on the color structures in which the final state SUSY 
particle pairs can be produced, contributes at the NLL 
level and is summarized by the factor 


Af\N, Q^,pA 


exp 


^ asCg) 

Jp q ^ 


( 6 ) 


The one-loop coefficients Dj follow from the thresh¬ 
old limit of the one-loop soft anomalous-dimension ma¬ 
trix l[35llJ7l. 

In Lig. [^we present the NLO-i-NLL predictions for 
squark and gluino production at the LHC {^^S = 8 
and 13TeV). The MS-scheme with five active flavors 
is used to define and the parton distribution func¬ 
tions (pdfs) at NLO. The masses of the squarks and 
gluinos are renormalized in the on-shell scheme, and 
the SUSY particles are decoupled from the running of 




Figure 1: NLO+NLL cross sections for squark and gluino production 
at the LHC (" VS =» and 13 TeV). (The FT curve in the lower plot 
corresponds to the t\t\ cross section; both notations q* and q are 
used in the literature to denote antis quarks.) The cross sections have 
been evaluated for a common squark and gluino mass mq = mg = 
^average = ^nd the renormalization and factorization scales have 
been set to p = m. The MSTW2008 pdfs ^67^ have been used. 


as and the pdf. The results have been obtained with the 
MSTW2008 parton distribution function IbVll . 

We observe that the inclusive cross section is domi¬ 
nated by qq and qg production at large sparticle masses, 
a consequence of the large valence quark component of 
the pdf at large x. Note that the qq, qg and qq cross sec¬ 
tions include the sum over five flavors of squarks, q e 
{fi, d, c, s, b}, and both chiralities {qp and qp). Therefore, 
and because of a different threshold behaviour, the cross 
section for the production of light stops, pp titi, is 
strongly suppressed. 

The NLL summation increases the cross section pre¬ 
diction if the renormalization and factorization scales 
are chosen near the average mass of the final state parti¬ 
cles. More importantly, threshold resummation leads to 
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Figure 2: The theoretical uncertainty for associated squark-gluino 
production at the LHC with 7TeV, pp qg + X, as a function of 
the sparticle mass mq = mg = m. The error bands represent the 
NLO+NLL scale uncertainty in the range mjl < p< 2m, and the total 
theory uncertainty including the 68% C.L. pdf and as error, added in 
quadrature, and the error from scale variation in the range mjl < 
p < 2m added linearly to the combined pdf and as uncertainty. The 
total theory uncertainty is shown at both NLO and NLO+NLL. The 
MSTW2008 pdfs ^7^ have been used. 


a significant reduction of the scale dependence over the 
full range of sparticle masses, with an overall scale un¬ 
certainty at NLO-rNLL of less than 10%. This is shown 
in Fig.[^ for squark-gluino associated production. We 
also show the full theory uncertainty, consisting of the 
68% C.L. pdf and as error added in quadrature, com¬ 
bined linearly with the scale variation error for the NLO 
cross sections and for the NLO-hNLL cross sections. We 
find that even though the pdf uncertainty is significant, 
the inclusion of threshold resummation leads to a size¬ 
able reduction of the overall theory uncertainty. 

The NLO-rNLL predictions for degenerate squark 
masses can be computed with the program Nll-fast |[68ll 
and are employed by the LHC experiments to interpret 
search limits and constrain the MSSM parameter space. 
Selected results for squark and gluino cross sections 
at current and future hadron colliders are collected in 
Refs, ismsa, including a more detailed discussion of 
the theoretical scale and pdf uncertainty. 

2.2. NNLL threshold resummation 

Let us first briefiy comment on the ingredients needed 
to extend the resummation of threshold logarithms to 


next-to-next-to-leading logarithmic accuracy. As dis¬ 
cussed in Section |2.1| the (soft-)collinear radiation ef¬ 


fects are described by the functions A/ and A^, while 
wide-angle soft radiation is accounted for by ^ in 
Eq. 0.. The product of these radiative factors can be 
written schematically as 

a,a,a‘®.^ ,,, = 

^ J I j^klJ 

exp [Lgi(asL) + giia^L) + asgsCasL) + (7) 

This exponent contains all the dependence on large log¬ 
arithms L = In N. The leading logarithmic approxima¬ 
tion (LL) is represented by the gi term alone, whereas 
the NLL approximation requires additionally including 
the g 2 term. Similarly, the g 3 term is needed for the 
NNLL approximation. The expressions for the gi and 
g 2 functions can be found in e.g. fSSj and that for the 
NNLL g 3 function in e.g. 14^ . 

The matching coefficients Cij^kij in Eq. ^ collect 
non-logarithmic terms as well as logarithmic terms of 
non-soft origin in the Mellin moments of the higher- 
order contributions. The coefficients Cij^kij factorise 
into a part that contains the Coulomb corrections and a 
part containing hard contributions | 




ij^klj 


(1 + -c 


s^CouUl) , ®s^Coul,(2) 


i j^klj 


^ j^2'^iHkl,I 


+ ...) 


a. 


(1 + -C'!' „,+ %C4>,,, + , 

^ ^ ij^klJ ^2 ij^klJ 




( 8 ) 


Apart from the terms of 0(as), which need to be in¬ 
cluded in Cij^kij when performing resummation at 
NNLL, some of the 0(a^) terms are also known and 
can be included in the numerical calculations. Expand¬ 
ing Eq. ^ we have 


C 


-NNLL _ 
ij^kU - 

l + ^(cCouUi) + c<i) 

JJ. V Ij^klj ij^kl, 


j 


^ /pCoul,(2) p(2) p(l) pCoul,(l)\ 

_2 y^ij^klj ^iHkl,I ^iHkl,I^iHkl,I ) • 


(9) 


The first-order hard matching coefficients ^ were 
calculated in 1491 , whereas the expressions for the first- 
order Coulomb corrections in Mellin-moment 

I j^kl,I 

space are listed in lEol. The form of the two-loop 
Coulomb corrections in x-space is known in the liter¬ 
ature 1691 . and the coefficients have been cal- 

culated in 1501 by taking Mellin moments of the near¬ 
threshold approximation of these two-loop Coulomb 
corrections. The second-order hard matching coeffi- 
cient C. j is not known at the moment and we put 

C®«,/ =0inEq. 









M. Kramer and M. Muhlleitner 


5 


Once we have the NNLL resummed cross section in 
Mellin-moment space, we match it to the approximated 
NNLO cross section, which is constructed by adding 
the near-threshold approximation of the NNLO correc¬ 
tion ||6^ to the full NLO result ifTTll . The matching is 
performed according to 


(NNLL matched)/ r 2^ 2\ (NNLOApprox)/ r 2i 2\ 


h\h2^kl 

5 i .= 




r ~ (res,NNLL) ^ 
ij^kl 


_ ~ (res,NNLL) ^ 


(NNLOApprox: 


,1. (10) 


To evaluate the inverse Mellin transform in Eq. ( p^ we 
again adopt the prescription of reference 1661 for the in¬ 
tegration contour CT. 

In Fig. we present the NNLL matched cross section 
prediction for the sum of the different squark and gluino 
production processes. The theoretical uncertainty in¬ 
cludes the scale error as well as pdf and errors. It 
is obtained by linearly adding the scale dependence in 
the range m/2 < fi < 2m to the combined 68% C.L. 
pdf and as uncertainties, the latter two added in quadra¬ 
ture. We see uncertainties grow from approximately 5% 
for masses near the present lower bounds, to 10% for 
masses approaching 2.5 TeV. Including the NNLL con¬ 
tributions leads to a further reduction of the scale depen¬ 
dence for the squark and gluino production total cross 
sections, with the exception of gluino-pair production 
as discussed in more detail in Ref. l50l . 


2.3. NLO-QCD corrections for generic MSSM spectra 
The calculations implemented in Prospino and Nll- 
FAST described in the previous section assume that five 
flavors of left- and right-chiral squarks, d^R, Cl,r, 
sl,r, and bL,R, are mass-degenerate. NLO-QCD predic¬ 
tions for generic MSSM spectra have been presented in 
the literature recently lESEll, and we shall briefly com¬ 
ment on these calculations below. 

We shall follow Ref. 1571 and discuss the production 
of squark-antisquark pairs as an example. General NLO 
squark and gluino production cross sections have also 
been presented in |[54l . At LO squark-antisquark pro¬ 
duction, pp qq, proceeds through quark-antiquark 
annihilation and gluon-fusion, 

~cl ~c2 

’ (11) 

88 • 

Here, the lower indices indicate the flavor of the par¬ 
ticle, whereas the upper indices for the squarks denote 


10^ 
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Figure 3: The NNLL matched cross section for the sum of the four 
processes of pair production of squarks and gluinos, including the 
theoretical error. The error band includes the 68% C.L. pdf and cTs 
errors, added quadratically, and the scale uncertainty varied in the 
range may 12 < p < 2mav, added linearly to the combined pdf and cTs 
error. The energy is that of the LHC at 8 TeV. The squark and gluino 
masses have been taken equal and the common renormalization and 
factorization scale has been set equal to the average mass of the two 
particles produced. The MSTW 2008 NNLO parton distribution func¬ 
tion Il67\ has been used. 



the chirality. Due to the flavor conserving structure of 
the relevant vertices the gg initiated processes and the 5*- 
channel gluon exchange only contribute to the produc¬ 
tion of squarks of the same flavor and chirality. We shall 
focus on the production of squarks of the first two gen¬ 
erations. In total, this includes 64 possible final state 
combinations. This number can be reduced to 36 inde¬ 
pendent channels if the invariance under charge conju¬ 
gation is taken into account. 

To illustrate the size of the SUSY-QCD corrections 
for the individual production channels, let us con¬ 
sider a particular benchmark scenario of the constrained 
MSSM fTOl . specified by universal GUT-scale param¬ 
eters mo/mi/ 2 Mo = 825/550/0 GeV, tan(yS) = 10 and 
sgn(p) = -rl. The corresponding squark and gluino 
masses read 


muL = ^CL 

II 

II 

II 

m~g 

1799.53 

1760.21 

1801.08 

1756.40 

1602.96 


where we have neglected the very small mass splitting 
between the first and second generation squarks. 

In Table we present LO and NLO results for in¬ 
dividual squark-antisquark production channels at the 
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Process 

CLoLfb] 

c^'nloM 

K-factor 

UlUl 

9.51 

10-2 

1.43 

10-1 

1.50 

UrUr 

1.14 

10-1 

1.72 

10-1 

1.51 

didi 

5.50 

10-2 

8.79 

10-2 

1.60 

dRdR 

6.89 

10-2 

1.11 

10-1 

1.61 

UlUr 

3.75 

10-1 

5.12 

10-1 

1.37 

didR 

1.41 

10-1 

1.70 

10-1 

1.21 

uidi 

6.98 

10-2 

7.89 

10-2 

1.13 

UidR 

2.98 

10-1 

3.54 

10-1 

1.19 

Ur^l 

2.94 

10-1 

3.49 

10-1 

1.19 

URdR 

8.36 

10-2 

9.54 

10-2 

1.14 

Sum 

1.59 

2.07 

1.30 


Table 1: The LO and NLO cross sections for squark-antisquark pro¬ 
duction of the first generation at the LHC with -y/S = 14 TeV obtained 
for a specific scenario of the constrained MSSM as described in the 
text. The renormalization and factorization scales have been set to the 
average squark mass, i.e. = 1779.31 GeV, and the CTEQ6L1 (tT] 
and CTIONLO (72) parton distribution functions have been adopted 
for the LO and NLO cross sections, respectively. 


LHC with 14 TeV. The renormalization and factoriza¬ 
tion scales have been set to the average squark mass, 
i.e. = 1779.31 GeV, and the CTEQ6L1 IHl and 
CTIONLO 117211 parton distribution functions have been 
adopted for the LO and NLO cross sections, respec¬ 
tively. 

Note, that the channels with squarks of the same fla¬ 
vor and chirality in the flnal state, displayed in the first 
four rows of the table, have contributions from gg initial 
states and therefore larger .^f-factors than channels with 
squarks of different flavor or chirality. Hence, approxi¬ 
mating the individual .^-factors by the average .^-factor 
as obtained from Prospino, K^yg = 1.39, is not very ac¬ 
curate . 

Determining the individual corrections consistently 
is especially important if the decays are taken into ac¬ 
count and the branching ratios of the different squarks 
differ significantly for the specific decay channel under 
consideration. In order to assess the possible numeri¬ 
cal impact of this approximation we consider the decay 
q qx\ at LO at the level of total cross sections, i.e. we 
multiply the production cross sections for the individ¬ 
ual squark-antisquark channels with the respective LO 
branching ratios: 

^ crNLO-BRLO(§^;^o^).BRLO(§^;e??) 

channels 

( 12 ) 


Multiplying the LO result for each subchannel with 
the average ^-factor, K^yg = 1.39, and the correspond¬ 
ing branching ratios gives 

E ctlo • ■ BR^-O {q ^ ■ BR^o (§ ^ 

channels 

= 0.126fb. (13) 

Thus the rate obtained with the approximation relying 
on an average K-fdiCiox for all subchannels is roughly 
10% smaller for this particular case. 

To summarize this section, we have seen that using 
an average .^-factor to globally rescale individual sub¬ 
channels for squark-antisquark production may not be 
a very accurate approximation. For a generic bench¬ 
mark scenario, we And differences of (9(10%), i.e. of 
phenomenological relevance for the experimental accu¬ 
racy expected in run 2 of the LHC. In general, the more 
the branching ratios of the different squarks differ, the 
more important the consistent treatment of the individ¬ 
ual corrections becomes. 

3. Light stop decays 

Higher-order QCD corrections to supersymmetric 
particle decay widths and branching ratios have been 
calculated for many years, see e.g. Ref. |[73||74l and ref¬ 
erences therein. Here, we shall focus on the decays of 
light stops, which play a special role in supersymmetric 
models. 

In most SUSY models a light stop arises naturally, 
as the mixing is proportional to the large top Yukawa 
coupling, which leads to a large mass splitting between 
the stop mass eigenstates. Light stops play an impor¬ 
tant role in view of the Higgs mass and naturalness ar¬ 
guments; together with the top loops they provide the 
dominant higher-order corrections to the light CP-even 
Higgs mass ll75ti79]| , which are necessary to shift the 
tree-level mass to the observed value of about 125 GeV. 
A light stop can also lead to the correct dark matter relic 
density through co-annihilation, in particular for mass 
differences between the stop and the lightest neutralino 
x\ of 15-30 GeV |[80l - [85]| . Moreover, light stops allow 
for successful baryogenesis within the MSSM 1861 - 1981 . 

While the canonical LHC searches for jets and miss¬ 
ing transverse energy place lower limits on the masses 
of the squarks of the first two generations of approx¬ 
imately 1.5 TeV [131, the lightest stop can still 
be rather light and have a mass below the kinematical 
threshold for the decay into a top quark and the light¬ 
est neutralino If we assume the lightest stop 0 to 
be the next-to-lightest SUSY particle (NLSP) and the 


= 0.139fb. 
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lightest neutralino to be the lightest SUSY particle, then 
h can decay into the LSP and a charm quark c or an up 
quark u, i.e. t\ {ulc)x\ 1991IIOOII . Another possi¬ 
ble decay channel is the four-body decay ti x\^ff 
(Ml, where / and /' denote generic light fermions. 
The two-body decay into charm/up and neutralino is 
flavor-violating (FV). While the MSSM in general ex¬ 
hibits many sources of flavor violation, so that the decay 
can already occur at tree-level, high precision tests in 
the sector of quark flavor violation and limits on flavor- 
violating neutral currents from K, D and B meson stud¬ 
ies put stringent constraints on the amount of possible 
flavor violation. The hypothesis of Minimal Flavor Vi¬ 
olation (MFV) 11102111061 has been proposed in order to 
solve this New Physics Flavor Puzzle. It requires that ah 
sources of flavor and CP-violation are given by the SM 
structure of the Yukawa couplings, so that flavor mixing 
in models of New Physics is then always proportional 
to the off-diagonal elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix 111071 llOSll . However, the hy¬ 
pothesis of MFV is not renormalization group invariant 
ifTost . Flavor off-diagonal squark mass terms are in¬ 
duced through the Yukawa couplings, so that the squark 
and quark mass matrices cannot be diagonalized simul¬ 
taneously any more and e.g. the stop state receives some 
admixture from the charm (up) squark, inducing a tree- 
level flavor-changing neutral current (FCNC) coupling 
between the stop, the lightest neutralino and the charm 
(up) quark. Note that for very small FV stop-neutralino- 
up/charm quark couplings, the four-body decay can be¬ 
come important and has to be considered for a reliable 
prediction of the t\ branching ratios. The stop masses 
have been bounded by LEP 11091111011 and Tevatron 
mill I112L and more recently by ATLAS BllSll and 
CMS mi4L with the strongest limits coming from the 
ATLAS analyses Refs. 111511116ll . Ah these analyses 
assume a branching ratio equal to one for the consid¬ 
ered decay channel of the 0, i.e. either the FV two-body 
or the four-body decay. The branching ratios, how¬ 
ever, can deviate significantly from one in large parts 
of the allowed parameter range, as has been shown in 
Refs. mOOl 111711 and will be outlined in the following. 
Taking this effect into account, the experimental exclu¬ 
sion limits on the stop, which are based on the assump¬ 
tion of branching ratios equal to one, are considerably 
weakened. 

In Ref. moon we improved upon the existing approx¬ 
imate result for the t\ (u/c)x^ decay of Ref. ll99]| . 
which only takes into account the leading logarithms of 
the MFV scale. We calculated the exact one-loop decay 
width in the framework of MFV, by performing the full 
renormalization program and including the finite non- 


ii cx^i 

ri-ioop[GeV] 

r”/^[GeV] 

Scenario (1) 

9.322 • lO-'o 

1.004- 10-® 

Scenario (2) 

5.862 • 10-® 

6.446 • 10-® 


Table 2: The partial widths for the decay h cx^ in two MFV 
scenarios, calculated with the exact 1-loop formula, and with 

the approximate formula of Ref. (99l, F^/^. 

logarithmic terms arising from the loop integrals. For 
the numerical analysis we investigated two constrained 
MSSM scenarios with universal soft-breaking terms at 
the GUT scale Mqut ~ 10^^ GeV. This scale is iden¬ 
tified with the MFV scale, and ah soft SUSY breaking 
parameters are family universal. The boundary condi¬ 
tions at /iMFV = Mqxt are 

(1) mo = 200 GeV; mi /2 = 230 GeV; 

Ao = -920 GeV; tmijS =10; sign(/i) = -f; 

(14) 

(2) mo = 200 GeV; mi /2 = 230 GeV; 

Ao = -895 GeV; tanyS = 10; sign(//) = -f. 

The mass spectra and mixing angles have been calcu¬ 
lated with the spectrum calculator SPheno 11118111191 , 
which allows for flavor violation. The second scenario 
has a larger ii - x^^ mass difference compared to sce¬ 
nario (1), whereas the mass difference between ii and 
the lightest chargino | is smaller: 

(1) = 104/92/175 GeV; 

(2) = 130/92/175 GeV. 

The partial stop decay width into charm and neutralino 
has been calculated with the full one-loop formula and 
compared to the approximate result ll99ll . For the lat¬ 
ter, the charged W boson mass Mw has been taken as 
generic loop particle mass. The widths are given in Ta¬ 
ble]^ They have been obtained with the program Susy- 
Hit 117311741 . where the full one-loop formula for the 
flavor changing stop decay has been implemented. The 
results for the exact and the approximate decay widths 
differ by 0(10%). In the calculation of the branching ra¬ 
tios the ti ux\ and the four-body decay widths have 
also been included in the total width. The branching ra¬ 
tios are summarised in Table[^ In scenario (2) the larger 
mass difference between ti and x\ leads to a more im¬ 
portant four-body decay width, as it is dominated by the 
chargino exchange diagrams IIIQIII . This in turn induces 
a deviation of the branching ratio BR(fi ^ x^^c) from 
one by a few per-cent, as can be inferred from Table 

For MFV scales large compared to the EWSB scale, 
the two-body decay width is dominated by large log¬ 
arithms of the ratios of these scales, which should be 
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Figure 4: Exclusion limits in the m~o - plane at 95% CL, based on the results for the t\ cx\ signature from 11151 (left) and on the results 

for the t\ X^bf f' signature from I115II116I (right). The color code indicates the branching ratio down to which the exclusion limits are valid. 
Taken from 11171 . 


BR 

BR(fi ^f^c) 

BR(Fi 

Scenario (1) 

0.9944 

4.587 • 10-5 

Scenario (2) 

0.9443 

0.0504 


Table 3: The ii branching ratios for different final states for scenario 
(1) and (2). For details, see (Tool. 


resummed. Resummation effects can be included by 
taking into account the FCNC couplings that are in¬ 
duced through renormalization group running already 
at tree-level. This has been done in Ref. II117II where we 
calculated the one-loop SUSY-QCD corrections to the 
two-body decay. For the correct determination of the 
ti branching ratios, also the four-body decay is com¬ 
puted by consistently including FCNC couplings. Fur¬ 
thermore, non-vanishing masses for the third genera¬ 
tion fermions in the final state have been taken into ac¬ 
count. These decay widths have been implemented in 
Susy-Hit fTSllT?! . (The program with the newly imple¬ 
mented stop decays is available at lEol .) The improved 
ii branching ratios have implications for the LHC stop 
searches and the bounds obtained on the mass of the 
lightest stop . To show this we have performed a pa¬ 
rameter scan and taken care to respect constraints from 
Higgs data, SUSY searches, 5-physics measurements 
and from the relic density. For details, we refer to II117II . 
Figure shows the exclusion limits on depending 
on the mass m^o, taking into account our results, i.e. as 
a function of the actual value of the ti branching ra¬ 
tio. The plot reinterprets the ATLAS searches based 
on monojet-like cia and charm-tagged event selec¬ 


tions nisi and on searches for final states with one 
isolated lepton, jets and missing transverse momentum 
ESI, where limits on the lightest stop mass as a func¬ 
tion of the neutralino mass have been given, assuming a 
branching ratio of one. The grey dashed lines in the 
figure limit the region in which m^o -f < ruf^ < 
+ mw • In this region the stop can be searched 
for in the FV two-body decay and the four-body decay. 
The full pink line in the upper plot is the 95% CL exclu¬ 
sion limit based on combined charm-tagged and mono¬ 
jet ATLAS searches in the ti cx\ decay I115L as¬ 
suming a branching ratio of one. Under the assumption 
that ?i decays exclusively into the four-body final state, 
ATLAS derived from the monojet analysis ES the ex¬ 
clusion given by the pink line (close to the upper dashed 
line) in Fig. (right) and from the final states with one 
isolated lepton the exclusion region delineated by the 
green line (close to the lower dashed line) ifTTbl . With 
the information given in 1115111161 we derived the ex¬ 
clusion limits for the two- and the four-body final state 
as a function of the branching ratio, which is given by 
the color code. 

The plot shows the stop masses that are excluded for 
a branching ratio above the one associated with a spe¬ 
cific color. Evidently, for smaller branching ratios the 
exclusion limits become weaker. The combination of 
the two plots allows to extract the exclusion limits for 
stops of a given mass as function of the neutralino mass 
and the stop branching ratio. Thus it can be read off 
from Fig. (left) that ?i masses of 150 GeV can be 
excluded for m^o = 80 GeV if their branching ratio 
into c + x\ exceeds 0.43. This in turn implies that the 
stop four-body branching ratio is below 0.57. On the 
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Other hand the right plot shows that in the same region 
stops can be excluded if their branching ratio into the 
four-body final state is larger than 0.88, which implies 
that the two-body decay branching ratio is below 0.12 
then. This means that =150 GeV can be excluded 
for m^o = 80 GeV for scenarios in which BR(fi ^ 

A 1 

< 0.12 and BR(?i ^ cx^) > 0.43, respectively, 
BR(fi -^Xibff) > 0.88 andBR(fi ^ f,bff) < 0.57. 
The dark blue regions correspond to stops with vanish¬ 
ing branching ratios, so that all stop mass values asso¬ 
ciated with these regions are excluded. In Fig. (left) 
there is no smooth transition between the dark blue and 
its neighbouring regions, as the exclusion limits in the 
two-body final state are related to the ones in the four- 
body final state which here apply for branching ratios 
> 0.46, so that of course also in Fig. (right) there is 
no continuous color gradient. 

The exclusion limits given by the border of the col¬ 
ored region at 100% two-, respectively, four-body decay 
branching ratio, do not exactly match the ones derived 
by ATLAS. The reason is that ATLAS provided infor¬ 
mation on the values of the excluded production cross 
section times branching ratio only for a few points in 
the - mfj plane so that a linear interpolation be¬ 
tween these points was necessary in order to cover the 
whole region. Nevertheless, the agreement of the pre¬ 
sented results with the given exclusion limits is reason¬ 
ably good. The advantage of our approach is that it 
takes properly into account the information on the actual 
stop branching ratios which can considerably weaken 
the stop exclusion limits as is evident from Fig. As 
these plots can only be an approximation of what can 
be done much more accurately by the experiments, they 
should be taken as an encouragement to provide results 
also as function of the stop branching ratios. 

4. Exclusive squark production and decay at the 
LHC 

Fully differential NLO calculations of squark and 
gluino production and decay have been presented re¬ 
cently (531 I55H571 11211 I122L including the match¬ 
ing to parton showers for squark-squark and squark- 
antisquark production ESIEtI. These calculations show 
that higher-order QCD effects may significantly modify 
the shape of differential distributions. Thus, leading- 
order Monte-Carlo predictions scaled with inclusive 
QCD corrections do, in general, not provide an accurate 
description of exclusive observables and cross sections 
with kinematic cuts. We shall discuss the impact of the 
higher-order corrections on the differential distributions 



Figure 5: Transverse momentum distribution of the hardest jet, , 
for squark production at the LHC (^|S = 14 TeV), combined with the 
subsequent (anti-)squark decay qjq qlqx\ for the SUSY scenario 
described in the text. Shown are the LO and NLO results, as well as 
the differential ^-factor (full) and the total ^-factor (dashed). 

following the analysis of squark-antisquark and squark- 
squark production presented in Ref. Ell. 

In order to quantify the impact of the parton shower 
on the differential distributions, we have matched the 
NLO-calculation presented in Section |2.3| with various 
parton shower showers using the matching scheme of 
Ref. II 12311 as implemented in the Powheg-Box II124L 
The corresponding code is publicly available and can 
be downloaded at (na. 

Let us first comment on the impact of the NLO cor¬ 
rections on the shape of the differential distributions. 
We consider the production of a squark-antisquark pair, 
pp qq, followed by (anti-)squark decay into a neu- 
tralino and a jet, qjq qlqXy Fig. ^ we present 
the LO and the NLO distributions for the transverse mo¬ 
mentum of the hardest jet, p ^, using the benchmark sce¬ 
nario described in Section [2?^ (see Ref. (5711 for details). 
One observes a strong enhancement of the NLO correc¬ 
tions for small values of pr, while they turn even neg¬ 
ative for large values. The result for the second hardest 
jet, which is not shown here, is qualitatively the same. 

More accurate predictions for differential distribu¬ 
tions and exclusive observables require the match¬ 
ing of the NLO calculations with a parton shower. 
We have followed the matching scheme presented in 
Ref. II123II and implemented the NLO SUSY-QCD cal¬ 
culation of pp (q ^i)(§ ^ W?) + ^ iiito 
the Powheg-Box II124L We have compared the default 
Pythia 6 lfT26l and Herwig-i--i- lfT27l parton showers, 
and a pr-ordered dipole shower 1112811129L which is 
also implemented in Herwig-i-i-. We find that the pre¬ 
dictions of the different parton showers for the observ¬ 
ables depending solely on the two hardest jets agree 
within 0(10%) or better. Comparing the showered re- 
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Figure 6: Transverse momentum distribution of the hardest (upper 
panel) and third hardest (lower panel) jet, for squark production at the 
LHC (^/S = 14TeV), combined with the subsequent (anti-)squark 
decay qlq qlqx\ for Ihe SUSY scenario described in the text. The 
NLO predictions and the NLO results matched to the parton showers 
Pythia Q 26 l, Herwig++ and the Dipole-Shower mu [EH, 

are shown. The lower inserts depict the respective ratios of the results 
obtained with the three parton showers and the pure NLO prediction. 

suits with the outcome of a pure NLO simulation the 
effects of the parton showers on these observables are at 
most of (9(20%), except for the threshold region. Larger 
deviations between the different parton showers emerge 
in the predictions for the third hardest jet, which is for¬ 
mally described only at LO in the hard process, see 
Fig.0 

Let us finally comment on the impact of the higher- 
order corrections on event rates after experimental se¬ 
lection cuts. We consider the definition of a signal re¬ 
gion for the SUSY searches in two-jet events performed 
by the ATLAS collaboration with event selection 
cuts corresponding to 

> 130 GeV, > 60GeV, Bt > 160 GeV, 

— > 0.2, > 1 TeV, Ir) > 0.4, 

Weff 

and A^O- 3 ,lr)> 0.4 if p^^>40GeV. (16) 

Here, the effective mass meff is defined as the sum of the 
Pt of the two hardest jets and f j, whereas the inclusive 



qq 


NLO 

0.871 fb 

0.0781 fb 

Pythia 

0.883 fb 

0.0797 fb 

Herwig-i-i- 

0.895 fb 

0.0807 fb 


Table 4: Total cross sections after applying the event selection cuts 
defined in Eq. d for squark-squark and squark-antisquark produc¬ 
tion and subsequent decay q qx\ within the benchmark scenario 
introduced in Section [23] and specified in detail in Ref. ED Results 
have been obtained at the level of a pure NLO simulation and includ¬ 
ing parton shower effects with Pythia and Herwig-i--i-, respectively. 

definition of this observable includes all jets with p:! > 
40 GeV, 

Ki = ^4+BT. ( 17 ) 

i=l 

Moreover, A 0 ( 7 /,fj) denotes the minimal azimuthal 
separation between the direction of the missing trans¬ 
verse energy, tr, and the jet. The additional cut 
> 0.4 is only applied if a third jet with 
p^r^ > 40 GeV is present. 

Applying these cuts at the level of a pure NLO simu¬ 
lation yields the event rates as shown in the first row of 
Tab.|^ Results are shown for squark-squark and squark- 
antisquark production and subsequent decay q qx\ 
within the benchmark scenario introduced in Section [2^ 
and specified in detail in Ref. liSTl . Matching these NLO 
results with a parton shower hardly affects the outcome 
after using the cuts defined in Eq. ( p^ , as can be in¬ 
ferred from the results obtained with Pythia and the 
HERWiG-f-h default shower listed in the second and third 
row, respectively. 

To conclude, we have discussed the effect of match¬ 
ing various parton showers to a NLO-QCD calculation 
of squark-squark and squark-antisquark production. We 
find moderate shower effects of at most of (9(10 - 20 %) 
for the observables depending solely on the two hardest 
jets, or for rather inclusive observables like event rates 
with selection cuts. In more exclusive distributions, in¬ 
cluding for example those for subleading jets, the effect 
of the parton shower may be more sizeable. 

5. Conclusions 

Supersymmetric theories are among the most attrac¬ 
tive extensions of the Standard Model, and searches for 
SUSY play a key role in the LHC physics program. 
A tremendous amount of work has been done in the 
last 20 years to calculate higher-order corrections to 
the production and decay of supersymmetric particles 
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at hadron colliders. Such precision calculations are im¬ 
portant to place limits on the masses of supersymmet¬ 
ric particles from current and future LHC searches, and 
they will be essential to determine the parameters of a 
supersymmetric theory in the case of discovery at the 
upcoming LHC run 2. 

The production cross sections for the strongly inter¬ 
acting SUSY particles, squarks and gluinos, are known 
at next-to-leading order in SUSY-QCD, including the 
summation of threshold logarithms at next-to-next-to- 
leading accuracy. The uncertainty of the theoretical pre¬ 
dictions is at a level of < 15% and dominated by the 
error on the parton distribution functions. Similarly, 
there is a vast literature on higher-order corrections to 
SUSY particle decay widths and branching ratios, cov¬ 
ering most of the relevant decay modes in the MSSM. 

To fully exploit the potential of the LHC run 2, it 
is crucial to obtain precise predictions for more exclu¬ 
sive observables, including differential distributions and 
cross sections with cuts. To this end, one needs to con¬ 
sider higher-order corrections to SUSY particle produc¬ 
tion and decay chains, and to match the NLO predic¬ 
tions with parton showers. First results of such a more 
comprehensive analysis have been presented here. 

Finally, future work should focus on developing auto¬ 
mated tools to extend precision calculations to a wider 
class of beyond the SM physics scenarios. Such tools 
need to combine next-to-leading order accuracy in QCD 
for generic beyond the SM processes with the summa¬ 
tion of large logarithmic threshold corrections and the 
matching with parton shower Monte Carlo programs. 
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